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A criterion for effective irrelevancy of the spin-orbit coupling in the heavy- fermion superconduc- 
tivity is discussed on the basis of the impurity Anderson model with two sets of Kramers doublets. 
Using Wilson's numerical renormalization-group method, we demonstrate a formation of the quasi- 
particle as well as the renormalization of the rotational symmetry-breaking interaction in the lower 
Kramers doublet (quasispin) space. A comparison with the quasispin conserving interaction exhibits 
the effective irrelevancy of the symmetry-breaking interaction for the splitting of two doublets A 
larger than the characteristic energy of the local spin fluctuation Tk. The formula for the ratio of 
two interactions is also determined. 
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There has been a long standing issue on strength of 
spin-orbit coupling since Anderson pointed out its impor- 
tance in heavy- fermion superconductivity 0|. In heavy- 
fermion systems, the bare spin-orbit coupling is con- 
siderably large compared to a renormalized bandwidth 
SB: ill which heavy- mass quasiparticle forms 
Cooper pairs0. Experimentally however, there are sev- 
eral contradictions to this naive expectation. In the 
triplet superconductor, UPta, the so-called d- vector can 
rotate freely in application to relatively low magnetic 
fields which suggests a weak pinning force due to 

the spin-orbit coupling. Moreover, the potential candi- 
dates of the triplet superconductors, such as UNi2 AI3 and 
URu2Si2, exhibit power-law temperature dependence in- 
dicating the presence of line of zeros 0, which is incon- 
sistent with the group theoretical argument based on the 
strong spin-orbit coupling [l|, |3, llO? lul • 

The key to resolve the inconsistency is a renormal- 
ization of the spin-orbit coupling inherent from a for- 
mation of the quasiparticle. Miyake phenomenologi- 
cally argued this point taking the lowest Kramers dou- 
blet in crystalline-electric-field (CEF) as a Wannier basis 
of the quasiparticle tight-binding model 12], where the 
bare strong spin-orbit coupling has been strictly taken 
into account. Provided that the primary quasiparticle 
band is well isolated from any other bands, residual rota- 
tional symmetry-breaking interactions are of the order of 
(£^p/A)^, A being the splitting of the bands. It ensures 
an approximate conservation of the rotational symme- 
try in the quasispin (Kramers doublet) space. Although 
the scenario seems to be plausible, it is highly non-trivial 
whether the simple one-body picture properly works for 
the quasiparticle and the CEF excited states or not, as 
we recall many-body aspect of the problem. 

In this paper, we demonstrate the renormalization of 
the spin-orbit coupling as well as the quasiparticle forma- 
tion, and elucidate the criterion for weak residual spin- 
orbit interaction in the heavy- fermion systems. For this 
purpose, we solve the impurity Anderson model with 



two set of Kramers doublets usi ng W ilson's numerical 
renormalization-group method [l^. UA . [l5 L Vwi [l7| . which 
is known to be numerically exact for impurity Ander- 
son models. Viewing a renormalized local interaction as 
a source of the momentum-dependent effective interac- 
tion for superconductivity [23, a comparison in strength 
among local interactions gives a useful threshold for an 
irrelevancy of the residual spin-orbit interaction. We as- 
sume that the crystal has the inversion symmetry, other- 
wise the quasispin degeneracy of the band will be lifted. 
In this case the spin-orbit coupling competes with the 
condensation energy rather than the attractive interac- 
tion leading to the superconduc ting state. It is a separate 
issue from the present paper [H llfj En|. 

We begin with the impurity Anderson model with two 
Kramers doublets. 
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Viflck^a + h.c.) 



(1) 

where we have neglected k and orbital (7) dependences 
in the hybridization and the dispersion of the conduction 
electrons. We restrict ourselves to the case of infinite 
Coulomb repulsions among all orbital and spin states for 
simplicity. As a simple example for two Kramers dou- 
blets, we adopt \T7±) = | ± 1/2) and iFgi) = | T 3/2) 
in the hexagonal symmetry (7 = 7 or 9). Since |m) 
represents the eigenstate of the 2;-component of the to- 
tal angular momentum j, the bare spin-orbit coupling is 
strictly taken into account at this stage. Then, the f^-f^ 
restricted /-electron Hamiltonian reads 



Ef\Tra){rra\ + {Ef + A)|r9a)(r9a| 



, (2) 



where A is the splitting of two doublets, and is the ad- 
justable parameter in the following discussion. We fix 
the other parameters as = — 1 and W = = 0.05 in 
unit of the half-width of the conduction band. 
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FIG. 1: The spectral intensity of the dynamical transverse 
susceptibility for several CEF splittings A at T = 0. For 
A > Tk (Tk is determined by the (lower) peak position), the 
second broad peak at cj ^ A is pronounced, otherwise the 
CEF excitation is obscured by the Kondo resonance scatter- 
ing. 



First, let us characterize (local) quasiparticle for sev- 
eral values of A. Figure shows the spectral intensity 
of the dynamical transverse susceptibility, Imx±(^) at 
T = 0. As A increases, the quasi-elastic peak becomes 
sharper with its position going downward. This is under- 
stood by the suppression of spin fluctuations between or- 
bitals, which gives rise to a reduction of the Kondo energy 
scale Tk- Here we have defined Tk at which lmx_L{uj) 
shows the maximum. When A exceeds Tk (in cases of 
A = 0.005, 0.01), the second broad peak at cc; ~ A corre- 
sponding to the CEF excitation is pronounced, otherwise 
the CEF excitation is obscured by the Kondo resonance 
scattering 2 Ij. 

The /-electron density of states for each orbitals 
pfj{uj) is shown in Fig.O Due to the particle- hole asym- 
metry in the model, the Kondo resonance (local quasi- 
particle) peak appears slightly above the Fermi energy. 
Reflecting the reduction of Tk with increase of A, the 
quasiparticle peak becomes sharper in the lower Tj or- 
bital, while the peak in the upper Fg orbital turns into 
the broad incoherent part of the spectrum located at 
CO ~ 2 A. Since with increasing A the integrated spec- 
tral weight below the Fermi level increases (decreases) in 
the Fy (Fg) orbital, the /-electron number changes to- 
ward unity (zero) (see also the inset of Fig. Note 
that the total /-electron number is almost independent 
of A, i.e., rif ^ 0.95. Therefore, we can conclude that 
the dominant component of the quasiparticle comes from 
the lower CEF doublet via the local spin fluctuation with 
conduction electrons, as long as the CEF excitation is 
clearly observed in the neutron scattering intensity. 

Next, we discuss the renormalization of the spin-orbit 
coupling. As was mentioned that the bare spin-orbit cou- 
pling has already been included in the definition of the 




FIG. 2: The /-electron density of states for Ft (left panel) and 
Fg (right panel) orbitals. With increasing A, the quasiparticle 
peak becomes sharper in the lower Fr orbital, while the peak 
in the upper Fg orbital turns into broad incoherent part of 
the spectrum. 

CEF state, we consider here two type of residual interac- 
tions, i.e., quasispin conserving and non-conserving inter- 
actions, respectively. To quantify the extent of the renor- 
malization, let us define the effective transition probabil- 
ity for an operator O as 

Po ^ ^e-'^(^'"+^")|(n|d|m)|V^', Z = ^e"'^^" 

mn n 

(3) 

where En and |n) are the exact eigen energy and its eigen- 
state, respectively, and /3 = 1/T is the inverse tempera- 
ture. Then, we consider Ointer = /g^/rt ^^e quasispin 
non-conserving operator, which is the essential part of 
the "spin-orbit" coupling, and Ointra = fj^fii foi" the 
conserving operator. Note that Po with Ointra (we de- 
note it as Pintra hereafter) should be equal to 2Po with 
O = (/^^/7| — fj^fj^)/2 due to the Pauli principle. 

The temperature dependence of Pintra and Pinter for 
A = 0.003 is shown in Fig.|31 The thermal average of /- 
electron number in each orbitals is also shown in the inset 
of Fig. El As we expect, the orbital polarization in the /- 
electron occupancy is pronounced for T < A. For T > A 
two probabilities are almost T independent, and have 
the same magnitude. In the region Tk < T < A, Pinter 
begins to decrease rapidly while Pintra remains T inde- 
pendent. With further decrease of T, Pintra also starts to 
decrease with the same exponent as that of Pinter- Then, 
the ratio of two probabilities, r(T) = Pinter/Pintra is unity 
in the region, T > A, while r(T) <C 1 and T independent 
for T < Tk, provided that A > Tk- In other words, when 
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FIG. 3: Temperature dependence of the effective transition 
probability of the intra- and the inter-orbital spin-flip oper- 
ators for A = 0.003 (see text in detail). The ratio of two 
probabilities is unity in the region, T > A, while it is consid- 
erably small and T independent for T < Tk, provided that 
A > Tk. The inset shows the thermal average of /-electron 
number in each orbit als. 

the quasiparticle is formed below Tk, and if Tk is smaller 
than the first CEF excitation energy, the quasispin non- 
conserving interaction is negligible as compared with the 
quasispin conserving interaction. Consequently, the pair- 
ing interaction for the quasiparticle has the approximate 
rotational symmetry in the quasispin space. 

Now, we discuss more accurate criterion for the irrel- 
evancy of the rotational symmetry-breaking interaction 
as a function of Tk and A. In general, the transition 
probability near the fixed point is expressed as 

Po-{TlT:^f^^, (4) 

where Tq is the characteristic energy of the operator 
O and do is its scaling dimension in the scale trans- 
formation near the fixed point 0, [2^. Near the high- 
temperature fixed point, the /-component in the wave- 
function is localized. Then, both operators, Ointra and 
Ointer tiavc the scaling dimension d = 0, yielding the 
T independent behaviors at high temperatures. For 
T < A, Pinter scts in the regime of the low-temperature 
fixed point, while Pintra stays in the regime of the high- 
temperature fixed point. Since the / wavefunction of the 
Fg orbital is delocalized in the low-temperature regime, 
the contribution from the impurity site decreases expo- 
nentially in the scale transformation. Indeed, the lo- 
cal hopping such as Ointer has the Fermi-liquid scaling 
dimension, d = 1. Then, Pinter displays the relation 
In(Pinter) ^ 2 ln(T/T*) for T < T* - A. For T < Tk, 
Pintra Gutcrs the low-tcmpcraturc regime as well, showing 
the similar T dependence as Pinter with T* ~ Tk in this 
case. In this way we can understand the T dependences 
obtained in Fig. 13 

We plot the zero-temperature limit of the probabil- 
ity ratio ro = r(T ^ 0) as a function of A in Fig. ^ 
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FIG. 4: The ratio of the intra- and the inter-transition prob- 
abilities as a function of the CEF splitting (square). The 
A dependence of the ratio is well described by the formula, 
(TK/(TK + aA))2 with a = 2.5 (cross). 



(square). From the above discussion, ro ~ (^intra/^inter)^ 
at low temperatures. If we adopt Ti*tra ^ and 
^inter + <^A, Oi being of the order of unity, we have 

ro ~ [Pk/(Pk + aA)]^. Indeed, the ratio of two prob- 
abilities is well described by the formula with a = 2.5 
as shown in Fig. ^ (cross). It should be noted that Tk 
itself strongly depends on A, namely, the increase of A 
suppresses Tk- As a result, the ratio decreases faster 
than A~^. Recalling the exponential dependence in the 
critical temperature of the BCS formula, quasispin non- 
conserving interactions are practically irrelevant to the 
Cooper-pair formation in the quasiparticle band. 

It is interesting to note that a similar situation occurs 
in Sr2Ru04. In this case, the bare spin-orbit coupling 
A rather than the renormalized one appears in the sec- 
ond order, A^/T^p, because the primary 7-band for the 
superconductivity is isolated from the other a/3-bands 
due to the symmetry reason 23j. It can be said that 
in the heavy-fermion systems, the process of the heavy 
quasiparticle formation dynamically provides an isolated 
quasiparticle band similar to the 7-band in Sr2Ru04. 

So far, we have assumed the /-electron valency is less 
than one. In this case, a characteristic of the quasipar- 
ticle in the lattice system and a process of its forma- 
tion are well understood^, 0, H- On the other hand, 
we have no concrete understanding for a situation of 
more than two /-electron valency. In particular, the /^ 
configuration with a singlet CEF ground state, which 
is believed to be realized in some U-based conipounds 
and Pr-based filled skutteruditesjH IH HI HI 111 , 
is remarkable. This is because the / electrons them- 
selves can release their entropy without hybridizing the 
conduction electrons US 1^ which is necessary to 
form the Fermi liquid for systems with the /^ config- 
uration. Such subjects are currently under extensive 
investigations!^ IH 113 El 113, El 113 
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In summary, we have studied the formation of the 
quasiparticle as well as the renormalization of the resid- 
ual spin-orbit interaction based on the impurity An- 
derson model with two set of Kramers doublets. Us- 
ing the numerical renormalization-group calculation, we 
have shown the effective irrelevancy of the rotational 
symmetry-breaking interaction in the quasispin space for 
the CEF splitting A larger than the local spin-fluctuation 
energy Tk- As the temperature decreases, the delocal- 
ization of the / electron first occurs in the upper CEF 
state. Subsequently, the / electron in the lower CEF 
state delocalizes to form the heavy-mass quasiparticle. 
As a result, the inter-orbital interaction, which generates 
quasispin non- conserving interaction, is smaller than the 
intra-orbital quasispin conserving interaction by a fac- 
tor of (Tk/A)^. Since Tk has strong A dependence, 
(Tk/A)^ ^ 1 is easily realized even for small A. There- 
fore, in most cases the residual spin-orbit coupling is ir- 
relevant for the heavy- fermion superconductors. 
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